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I. SYMBOLS 

4 

Dls of flight mechanics according to 
used: 

coefficients of the frequency equations 
Eurwitz determinants 

damping of oscillations in roll and yaw 

frequency of the oscillations in roll and 
yaw 

Froude number 
altitude 

control-bearing ratio 
measuring ratio (cf. p. 7) 

moment coefficients of the lateral 
stability (cf. pp. 12 and 13) 

aerodynamic time unit 
equilibrium flight Telocity 
variables of the frequency equation 

rotational angle about the vertical axis 
(yaw angle) * 

relative density of the mass of the airplane 

products of the pairs of roots of the 
frequency equation ' 

sums of the pairs of roots of the frequency 
equation 
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non— dimens ional time 


H. REASONS FOR AND AIM OF THE INVES TIGAIIGN 


For a longitudinally and laterally stable glide bomb design 
developed in the DFS "Ernst Udet" [l]^ two possibilities existed for 
obtaining a rectilinear flight path: either without automatic control, 

with maintenance of high manufacturing accuracy, in particular very 
good constructional symmetry, or else with the aid of an automatic 
control; in the latter case too great constructional accuracy could be 
foregone. Since, however, even slight constructional defects (probably 
unavoidable even in most careful manufacture) may lead to considerable 
disturbances [2], only the second possibility could be applied in 
practice. The longitudinal stability could be ensured very simply, 

(even without use of automatic control), and lateral stability, attain- 
able only with much more difficulty, became the main object; it was to 
be attained with the aid of a special control arrangement in the glide 
bomb. In order to limit to a minimum the drop tests which are expensive 
and require much time, the lateral— stability properties of the bomb and 
the most favorable design of the automatic pilot and the most favorable 
adjustment of the control possible within practical limits were to be 
clarified by calculation beforehand. This is, essentially, the aim of 
the present report. 


IH. DESCRIPTION OF THE AUTOMATICALLY— 
CONTROLLED GLIDE BOMB 


Figure 1 is an outline drawing of the investigated glide bomb. 

The fuselage consists of a bomb SC-^OO, the cylindrical part of which 
was slightly elongated; the short blunt tail section was replaced by a 
longer one of aerodyne, mically more favorable design, which was provided 
with a vertical fin and housed the control. The two wing halves, 
provided with ailerons, are attached laterally at the center of the 
fuselage. The wing shows pronounced sweep-bach but not dihedral and Is 
at zero incidence with respect to the fuselage. Since the wing section 
is symmetrical, both wing halves are of equal area. For attainment of 
a sufficient longitudinal stability and adherence to a certain c a ~ value 
an appropriate position of the center of gravity Is selected and, since 

^Numbers in brackets refer to the references at the end of the 
report. 
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tall surfaces are lacking, "both control surfaces are accordingly trimmed 
upwards. This then is the zero position about which they are differen- 
tially actuated to serve as ailerons. They may also he simultaneously 
deflected through a differential linkage as elevators . Chiefly for 
variation of the yawing moment due to sideslip (weathercock stability) 
a stabilizing fin was sometimes attached to the bomb nosej it slightly 
influenced, among other factors, the lateral force, the rolling moment 
due to sideslip, and the damping in yaw. 

The control required for attaining rectilinear flight and suffi- 
cient lateral stability was a gyro control developed in the DFS "Ernst 
TTdet" which operates without time lag, that is, effects any control- 
surface deflection pertaining to an angular displacement without delay. 
No device of any kind producing a leading of the control-surface 
deflection relative to the angular displacement was applied. Of course, 
any other control could be used provided it operated according to the 
same control law so that the results of this investigation apply to it, 
too. 


The control gyroscope is installed so that its measuring axis lies 
in the symmetry plane of the glide bomb and forms a tilt angle with the 
vertical axis; the purpose is sensitive measurement of banking and 
yawing deviations and their transformations into an aileron deflection 
through a gear with adjustable ratio. The proportionality between bank 
and yaw and aileron deflection is influenced by the gyroscope tilt 
angle. The following control method results: 

If, due to any disturbance, only a bank of the glide bomb occurs, 
a proportional aileron deflection counteracts it and after a few damped 
oscillations the "undisturbed— flight— path" position is re-established. 
An error occurring in the flight path? also is measured by the control; 
the originating aileron deflection produces a bank toward the side 
opposite to the course deviation. Thereby the flight path and with it 
(due to its weathercock stability) the glide bomb are, as intended, 
turned back in the direction of the undisturbed flight path (azimuth 
stability). With the approach toward the undisturbed flight path the 
bank is again neutralized so that an aperiodical transitional motion in 
the direction of the undisturbed flight path results which could, with 
reference to Mathias [ 3 ], be designated as "yawing" motion; damped 
rolling and yawing oscillations are superimposed on it. 


2 in case of calm air, this will apply throughout. 
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IV. SETTING DP OF TEE FREQUENCY" EQUATION 
(CHARACTERISTIC EQUATION) 
a. Preliminary Remarks 


For safety reasons a release velocity far "below the equilibrium 
velocity v 0 of the glide bomb was selected for releasing the glide 
bomb from the carrier plane so that in no case would an immediate 
ascent of the flight path occur. However, this leads to a path oscil- 
lation (phugoid) of the glide bomb; without suitable control measures 
it is improperly damped and has a relatively long period of oscillation 
so that during the considered time interval from the release of the 
glide bomb onward a steady flight path does not exist. Thus, investi- 
gations regarding the lateral stability of a glide bomb are of increased 
practical importance whan they take the simultaneous path oscillation 
into account; this will be attempted below. 

The numerical integration of the known six differential equations 
of the longitudinal and lateral motion, with consideration of the 
control law and with inclusion of terms of higher order, was not 
intended. Such a numerical integration would have required an intol- 
erable calculating expenditure and still would have given results in 
llmitedly accurate agreement with actual conditions, due to estimated 
aerodynamic coefficients of the glide bomb and to other simplifications . 
Accordingly a few approximations (suggested by the nature of the 
problem) were made which led to fan-reaching simplifications of the. 
calculation and to a no longer intolerable expenditure of work for the 
numerical evaluation. 

Of special Importance Is the assumption that c a , in this investi- 
gation practically dependent on a only, remains constant during the 
flight of the glide bomb. ThiB assumption is justified since the 
relatively properly damped ctr-o sc illation is considerably more rapid 
than the path oscillation, thus Is damped quickly, and a then retains, 
without change, approximately Its equilibrium value corresponding to 
the static longitudinal stability. Besides, this behavior — estimated 
by rough calculation — was well confirmed by film recordings of the 
releases. Hence all coefficients which are functions of a (or c a ) 
remain constant and lead to an extensive simplification of the 
calculation. 

Furthermore It Is presupposed, as is customary, that the course 
of the longitudinal motion is practically independent of the lateral 
motion. This Is justified here for the added reason that the coeffi- 
cients determining the longitudinal motion are, in the first approxi- 
mation, Independent of the unsymmetr leal 1 state of flight of the lateral 
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motion. This assumption is Justifiable only as long as one stays 
within the l imi ts of normal disturbances , the only ones considered here. 

The opposite, that the longitudinal motion has no significant 
effect on the lateral motion, no longer holds true here. This applies 
only when a generally steady and little disturbed flight path is being 
considered. Since, however, the deviations of the longitudinal motion 
from a steady flight path occurring in the present case far exceed the 
magnitude of small disturbances, they do influence the lateral motion 
quite noticeably and must be taken into account. In the complete 
equations of the lateral motion [4] there appear in connection with the 
longitudinal motioii, aside from the constant a, the values of v, 7, 
and the first derivative with respect to time, 7. It is true, these 
three values vary, but only with a frequency up to two orders of magni- 
tude lower than that of the rolling and yawing oscillation. 

If one, therefore, considers a small part of the course of the 
motion, sufficient for Judging the • stability of the lateral motion, 
after a disturbance, v, 7, and 7 vary so little in the meantime that 
they may be regarded as practically constant. 

This assumption enables solution of the three differential equa- 
tions of the lateral motion in the conventional manner without consid- 
erable difficulties. If one wants to find out whether lateral stability 
prevails during the entire course of the longitudinal motion, one sub- 
divides the course of motion appropriately and then introduces the 
corresponding values of v, 7, and 7 as sectionally invariable into 
the lateral motion. In the present investigation the lateral stability 
determined at the four points denoted in figure 2 was regarded as 
sufficient for the stability of the lateral motion for the full duration 
of a path oscillation and thus during the entire free flight of the. 
glide bomb. 

In order to Include with certainty even the largest path oscil- 
lations to be expected at the glide-bomb releases, the release velocity 
was assumed to be only 60 percent of the equilibrium velocity v Q of 
th'e glide bomb. 

The lateral stability for controlled longitudinal motion also may 
be investigated in the manner described above, although only under the 
limiting presuppositions that a remains practically constant or at 
least sectionally almo st constant and that the oscillation period of 
the path oscillation considerably exceeds that of the rolling and 
yawing oscillation. 

Furthermore the influence of the Mach number is neglected. Although 
a velocity of more than 200 meters per second, which corresponds to a Mach 
number of over 0.6, was attained in drop tests and though the present 
investigation is intended to include these velocities, a consideration 
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of- the Influence of the Mach number on the coefficients (which is known 
only very inadequately) would he so uncertain that we desisted from it. 
Besides, as tests hy Gothert [ 5 ] and [6] show, really critical condi- 
tions start appearing only at Mach numbers around 0.7; here, then, lies 
the limit of validity for the present investigation. 

Finally it should he mentioned that a change of the air density in 
the course of the lateral motion is not taken into consideration: 
first, the altitude of a single initial disturbance causing the oscil- 
lation (which is decisive for the determination of the lateral stability) 
does not vary so greatly as to render that profitable; second, the 
(after all, significant) difference in altitude between peak and bottom 
of the path oscillation (where the lateral stability is calculated) is 
included in the investigation of the influence of relative aircraft mass 
densities u B of different magnitude. 

Before we start setting up the equations for the lateral motion we 
shall discuss the control law of the built-in control since it is being 
introduced into these equations. 


b. The Control Law 

As described on page 4 , angles of rotation about the longitudinal 
as well as about the vertical axis are measured by means of the control 
used here, and a rolling moment proportional to them was produced 
through an aileron deflection without delay. However, since an aileron 
deflection generally also causes a yawing moment, this latter must be 
taken into account els well. The additionally possible slight variations 
of lift, drag, and lateral force are neglected since the aileron deflec- 
tion causing them remains small due to a disturbance assume dly small 
throughout. If the angle of rotation about the vertical axis to be 
designated as yaw angle is (for want of more suitable letters) denoted 
by k and is counted starting from the desired position, the control 
law reads: 


I = k( k + mqp) 


Therein k signifies the control— gearing ratio if one understands 
it as the ratio of the aileron deflection to the yaw angle with simul- 
taneous insensitivity to bank,* and m, the measuring ratio indicating 
what fraction of a bank is measured by the control like a yaw angle and 
transformed into an aileron deflection. 


^Translator’s note: Literally, "error— free bank." 
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In the arrangement selected here m depends on the degree of tilt 
of the control gyroscope in the glide tomb and remains throughout 
smaller than one, that is, the tank is always measured with less sensi- 
tivity than the yaw angle. 

Since later on, the rolling and yawing moments stemming from the 
control also are needed, the expressions for those are set up as well. 
Both moments are dependent only on ft and cp so that one may express 
them as follows : 

aileron rolling moment = L k *k + Ljp’cp 
aileron yawing moment = Nk*k + 

The moment essential for the lateral stability of the controlled 
glide tomb, newly added in comparison to the uncontrolled condition, 
is Iq,. It signifies an immediately effective restoring moment for the 
longitudinal axis which intervenes directly in case of tanking errors. 
This restoring ability about the longitudinal axis corresponds to the 
yawing moment due to sideslip which is important for the vertical axis 
and is regarded later in the same way as likewise characteristic vari- 
able. Taking the control law into account, one finds 1^’ = kmL|' 
which shows that not merely the control— gearing ratio, measuring ratio, 
or aileron effectiveness are of Importance, but their combined action. 

The moment L K * coming into effect in case of lateral deviations 
serves, as described on page 4 , only for maintaining the corn's e and 
thus is less inportant than I<p for the lateral stability of the glide 
bomb. 


The aileron yawing moment originating in case of an aileron actu- 
ation is produced only unintentionally. It does have a slight influence 
on the lateral stability as will be seen later; thus it must, as a 
precaution, be taken into consideration, but is otherwise of no 
importance. 

After these preliminary remarks the dynamic relations of the lateral 
motion are set up. The longitudinal and vertical axes of the glide bomb 
were selected as the two axes about which equilibrium of moments must 
prevail; for the equilibrium of forces a horizontal direction of refer- 
ence perpendicular to the flight path was chosen. The following sec- 
tion (c) offers the reasons for this latter choice which is at variance 
with convention. 
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o. Tile Horizontal Forces Perpendicular 
to the Flight Path 

Basically, any direction of reference vould be equally suitable 
for the equilibrium of forces if it only contains components of the 
forces responsible for the lateral path curvature. Preference is to be 
given to the one which appears especially illustrative from the dynam- 
ical viewpoint. 

According to what was said on page 6, the longitudinal inclina tion 
of the flight path y varies so little during the brief t ime interval 
sufficient for consideration of the lateral motion that it may be 
regarded as invariable. Thus the path section considered is part of a 
helix which one may visualize as wound on a vertical circular cy lind er 
of the radius of the flight path, as shown in figure 3. 

On this helical line travels the glide bomb which for the consid- 
eration with respect to forces of the lateral path equilibrium must be 
visualized as a mass point; the position of the body axis of the glide 
bomb is, for the time being, unimportant. The direction of the centrip- 
etal force which causes the lateral curvature of the flight path and 
coincides with the principal normal of the helix is horizontal and 
perpendicular to the flight path. Thus it is logical to balance the 
components of the air forces in this selected direction with the mass 
force which is effective as centrifugal force. All forces acting in 
the direction of flight path and gravity have no Influence on the 
lateral curvature of the flight path visualized as part of a helix, 
because they are perpendicular to its principal normal. 

Xn the specifications BUJ L 100 this axis used as the line of 
attack for the lateral-force equilibrium is not especially character- 
ized; however, it coincides with y g in case one puts x = 0. 

Actually (y„)_, „ is there perpendicular to z a as well as to x a 
ts X— 0 6 

and the cosines of the angles formed with the effective air forces 

may be taken directly from the attached form. Since the total air force 

is mostly given by its components A, W, and Q and since W, lying in 

the flight path x a , does not yield a contribution, the equilibrium 

ratio becomes: 


A cos(-z a ,y g ) x=0 + Q cos(y a ,y g ) Xss0 = J v(cos 7 )x 

o 


Additionally, it is here tacitly assumed that the laterally 
effective air forces — those stemming from the rotations of the glide 
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■bomb about its longitudinal and vertical axis as veil as from the 
aileron deflection — are so small that they may be neglected [3] and [7]. 

With A = c a jy v^F and Q = c^' £ v^FP and special consideration 

of the fact that Gr here must not be replaced by — £ y2j< ^ is 

cos 7 ^ 

customary in case of steady rectilinear flight, one then obtains 


c a £ v^F sin ji + c^’ ^ v 2 fP cos n = — 

S 


If wo / 


Since one has to deal for the lateral motion (however, not for the 
longitudinal motion!) only with small disturbances, one may 
equate sin u = u and cos p = 1 and there results 


°a f ^ ^ + c q* | ^0 = : T cos 7 x 
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( 1 ) 


d. The Moments about the Longitudinal Axis 

Here the complete Euler equation must be stated in order to mate 
consideration of the influence of the gyroscopic terms possible: 

L " J A + ( J z ~ Jy^z 


L contains all rolling moments stemming from the air forces; 
they are caused by sideslip, rolling and yawing, and by the aileron 
deflection. 

Since for glide bombs with relatively very small span the moments 
of inertia Jy and J z are almost equal, the slight difference between 

them, multiplied by the small rotational velocity cjuy stemming from 
the longitudinal motion and the yawing velocity, is so insignificant, 
especially in comparison with the strong aileron moment, that the right- 
hand gyroscopic term is neglected. Thus there remains after introduc- 
tion of al 1 partial moments 

V* 3 + V“s + L z ,co z + V<P + Lk * K = J A 

> , > < v ’ * tf ’ <1 y > 

sideslip rolling yawing control 
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e. The Moments about the Vertical Axis 


Euler's equation here reads 

E = J z O) z +(J y - 


Due to the very small span of the glide bomb J 2 is by about one 
order of magnitude smaller than Jy, and the expression — J x ^ may 

be replaced by J z . N is subdivided as in (d) and one obtains 


sideslip 


*— = — . — > + <- 
rolling yawing 


z’“z V < P +N b' K t /• \ 

wins? control ' ' 


(3) 


Here the gyroscopic term must not be neglectedj even though it is 
not Important , it is, compared to the yawing moments due to the air 
forces , sufficiently significant to be taken into consideration, in 
order to avoid errors. 


There is still another reason for the seemingly different evalu- 
ation of the gyroscopic terms in the moment equations for longitudinal 
and vertical axis : The sum of -the moments about the longitudinal axis 

is dominated by the aileron moment which acts on a large lever arm of 
about half the span and always may be made the predominant rolling 
moment by suitable design of the control and ailerons. Conditions are 
different for the yawing moments j there the air force yielding the 
restoring moment is not so widely variable as for the aileron, due to 
the spatial limitation of the vertical fin; furthermore, it acts on a 
lever arm which measures only centimeters and is, moreover, relatively 
uncertain due to manufacturing inaccuracies occurring in practice, 
caused either by the position of the center of gravity or other faults 
in constriction. Eor this reason the moments effective about the 
vertical axis (which is more sensitive) must be taken into consider- 
ation to a much higher degree than is required for the longitudinal 
axis, if surprises are to be avoided. 


f . The Solution of the Three Equations of Motion 
The equations of motion, here compiled once more 

c a § + Pq* £ v2l ’P = | r cos 7 x (l) 

Lp’P + + L z’ Q z + V<P + Lk ’ K = 


( 2 ) 
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Up* 3 + + V“z + iycp + N k ’k = J z (a> z + (3) 

are three simultaneous differential equations; however, their unknowns 
are not yet uniform. In itself, it would "be a matter of indifference 
which one is selected for the solution; however, since the first 
equation is particularly characteristic for the lateral motion, the 
three unknowns 4, 3 , and X appearing in that equation are to “be 
maintained throughout. 

Thus all variables appearing in the second and third equation must 
be expressed in terms of 4, 3 , and X and their derivatives with 
respect to time. This is done partly according to Eautenberg [8] and fa") 
with the presupposition that the disturbances of the lateral motion are 
small (thus the three unknowns just mentioned are only a little 
different from zero) and with the presupposition that for the longi- 
tudinal motion a and 7 — but no longer 7 — may be regarded as small 

quantities. Since steep nose dives may be left out of consideration 
anyhow, the following approximations, valid up to about 7 = -45°, 
result : 


9=4+3 tan 7 
K = 3 + X cos 7 


ajj. = 4 - x sin 7 
ay = 7 

o> z = 3 + X cos 7 


Furthermore, the moment increases L' and IT’, divided by the 
respective moments of inertia J x and J z , are expressed by the 
following non-dimensional moment coefficients of lateral stability: 


damping in roll 


rolling moment in yawing 


rolling moment due to sideslip 


aileron rolling moment (total) 
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aileron rolling moment (due to cp) 
abbreviatedly: aileron moment 



r_s i 2 8c L 

aileron rolling moment (due to 

*) 


= ‘ 

"s 1 2 Scl 

^xj S K 

rolling-yawing moment 



n 2 = 

"si 2 8cj$- V 
fzj s 

damp ing— in— yaw 



n z = - 

[si 2 dcir 2 

rzj Sn z £ 

yawing moment due to sideslip 
(weathercock stability) 


■ 

np =- 

[si 2 8c N 

LizJ 83 

aileron yawing moment (total) 



ng = 

■ s l 2 8c]j 

izj dfi 

aileron yawing moment (due to 

<P) 


> - [ 

J3_l 2 

H\ 8cp 

aileron yawing moment (due to 

i 

K) 


n K = j 

•si 2 Sc N 

TzJ 5 k~ 


For the four moment coefficients which are decisive for the 
restoring and damping moments about the longitudinal and vertical axis 
the sign was selected so that they appear — as generally custo mar y in 
the theory of oscillations — as positive quantities. 

Taking the control law valid here | = k(fc + mcp) into considera- 
tion, some moment coefficients may "be transformed, after "brief inter- 
mediary calculation, as follows: 

^ ' -“-e ~ A 

1|! = “ m n “ * ~ m if ^ 

By these transformations it was attained that the connection of 
the aileron moment 1^ which is particularly important for the lateral 
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stability with the known values of the control (k,m) and the automatic 
pilot (l|) becomes manifest and that the three remaining mome nt coeffi- 
cients are expressed by this aileron moment l m and a further known 

fn ft I 'r 

numerical value 


l1 i j 


2G- r T 

If, furthermore, the aerodynamic time unit Tc, = s = Dimension 

• covF L J 


and the relative aircraft mass density p = 

s 

results finally; 


20 - 

gpsF 


gpvF 

are introduced, there 


L*’ 


W 


Jx 

hL. 

Jx 


!•*' _ 



Dimension 



Dimension 


!/s 

Nx* _ 


l/s 

il z 
% z 

l/s 

2tL- 

J 2 

” 5* 

l/s 

Us 

T^ P 

l/s 2 

_ 

J z 

_ if£_ no 
T/ P 

l/s 2 

^S -1 

1/s 2 

w » 

V. = 

Jz 

_ ^ 1 

I| TjF ^ 

l/s 2 


l/s 2 

Nk 1 = 
Jz 

* 1| Tj- 2 ^ 

l/s 2 


Therewith the equations (l), (2), and ( 3 ), on pp. 11 and 12, 
arranged according to the three unknowns, p, p, and x* assume the 
following form: 


c a |i + c *3 - Tj, cos 7 X = 0 


(la) 
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1 , . 

+ ijy t 1 ? - (= + tanr )^] e + ^ M 

“ 5 !jT 2 M cos + ^C 3 * sin r + i 2 cos 7 )x 
+ (sin ?)X = 0 (2a) 

- + + tanr ) 1 d |3 - 4 - ? 

-;^iW° os * 

- - T^f) sin 7 + n z cos 7 jx -(cos ^X = 0 ( 3 a) 

F 

The further treatment of these three equations does not offer any 
additional peculiarities ; the unknowns, \i, 0, and X, are put propor- 
tional to e^ and X is then determined from the condition that the 
principal determinant of the previous equation system must vanish if 
the solutions are not to he identically zero. One then obtains for X 
an equation of the fifth degree, the characteristic equation: 


X^ + aX^ + bX-3 + cX^ + dX + e = 0 
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with the coefficients 

a = + n * + V " Ca tan 7 ) 

■b = H§_ i f + tan 7^ + 1 1 ^cp + n 0 

T F 2[Ll|\m / J 

+ i [lxn z - l z (sx - V) + (V" c * + n z)] 

c = ^({^[^ +taa ^ + X z] + ( n * “ ^)(5 + ta ^) 

*»,* + 1 ] ’■‘'P 

+ l^np - lp(ni - Tpr) + " n P tan r ) 

+ ^(c 4 '-c a tan7) [l x n 2 - l z (”x - ^)] J 

4 = i ‘ s (ii 1? + np X 

+' (c 1 >-c a te»^[^[^(| + te '‘ 7 ) + l «] 

+ ( %: “VX« + tan ^*“4^ 

+ o a |i p [n z ♦ - v)tonr]-(i, * 4 * tan ; 

•■»)>» 


e 
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Compared to the frequency equation of the uncontrolled lateral 
motion which is of the fourth degree, its degree is now increased by 
one. The reason lies in the fact that for the controlled glide bomb 
treated here, due to the control a formerly non— existent tie-up with 
a dir ection in space or, more accurately, with a directional plane, is 
added. 

X was not made non-dimensional by means of the aerodynamic time 
unit Tji because here that would not offer any advantage ; it would not 
make any difference for the later representation of stability domains, 
but would be of disadvantage for the determination of the oscillation 
frequencies since the result has to be converted again to seconds 
afterwards . 

Ho difficulty arises if, nevertheless, a non-dimensional representa- 
tion of the frequency equation should be required for some reason. One 
ha s to take into consideration that in solving for the three unknowns |i, 
P, and X one should have made them proportional, not to e^, but 
to e ZT , with the non-dimensional z appearing in the frequency 
equation instead of X and r in contrast to t, denoting a non- 
dimensional time. Since, therefore, zt must be equal to Xt and, after 
selection of an (at first arbitrary) time unit T, the non-dimensional 
time becomes ' T = i, there results X = k . ' If this value is inserted 
into the frequency 1 equation and the total equation is, moreover, enlarged 
by t5, it can be seen easily that due to the specific selection of the 
aerodynamic Tj. for T the Tj, powers in parentheses for all coeffi- 
cients cancel each other j one obtains the non-dimensional form of the 
frequency equation 

z 5 +. aTj-z^ + bTp 2 z3 + cTg.3 z 2 + dl-gAz + eTj.^ = 0 

It must be well considered here that — because of the nonsteady 
longitudinal motion due to a gyroscopic term — the Tji in the coeffi- 
cients bTp 2 , cTj-3 and dTy^ remains and thus in the non-dimensional 
frequency equation also the direct influence of the velocity is not 
cancelled. 


g. Brief Discussion of the Coefficients of the Frequency 
Equation (Characteristic Equation) 

As will be shown later (p. 20 ), one condition required for attain- 
ment of lateral stability is that all coefficients be positive. If a 
few numerical values, enumerated on pp. 25 and 26, are anticipated, one 
obtains for the separate coefficients under this point of view the 
following results: 
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a is with certainty always posit Its since the chiefly important 
dampings 1 2 and n z about the longitudinal and vertical axes are 
positive as is also the lateral— force increase and, in descending 
flight, — c a tan 7. Besides, this last expression equals c^. for 
steady rectilinear flight, but not under any other conditions. 

b also is positive because the restoring moments 1^ and no 
predominate for the longitudinal and vertical axisj due to the high 
relative aircraft mass density u s here — and also for the other 
coefficients — chiefly the terms containing that quantity are of 
importance. 

c depends on the expressions l^z and l x np and thus must 
always be made positive. 

d depends to a high 
of e. 

e will be transformed somewhat further; using the expression 
for on p.13 one obtains 

Similarly to the static lateral stability of uncontrolled flight 
this expression, too, has special significance. If one visualizes that 
the glide bomb sideslips to the side opposite to the bank during an 
error in banking, it is the rolling moment 1 <j, stemming from the 
aileron deflection and the weathercock stability m which attempts to 
right the glide bomb again whereas the rolling moment due to sideslip lp 
and the aileron yawing moment tend to increase the bomb's angular 

deviation. Thus sufficiently large values of and np must be 

selected in order to have e with certainty turn out positive, particu- 
larly if n® is known only approximately or the sign varies (in case 
of aileron deflections of different magnitudes), as happens occasionally. 

This brief discussion, which has shown that all coefficients of 
the frequency equation always may be made positive, will have to suffice 
for the time being; detailed treatment of the connections with the 
dynamic lateral stability will follow in later paragraphs. 


degree an 


■n 6 

Ig ^ + n P> 


thus on the sign 
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V. THE EVALUATION OF TEE FREQUENCY EQUATION 
(CHARACTERISTIC EQUATION) 
a. The Method Followed 


The frequency equation and knowledge of the values contained in 
it offer a means for making, in any case, a statement regarding the 
stability of the glide bomb after a disturbance of the lateral motion, 
and for calculating the occurring frequencies and dampings. This 
approach by particular values is practically the only one, if a 
comprehensive view of more general validity is to be given, as this 
investigation aims to do; a general solution is not possible and clear 
fundamental insight into the complicated inter-relationship of the 
variables of this problem is hardly ever obtainable. Thus it is the 
only possible course to strive for a comprehensive survey of the 
influence of the frequency equation on the lateral stability by numerical 
substitution of all variable quantities in automatic pilot and control 
and by the evaluation of the frequency equation. 

In limiting this investigation to an automatic pilot of the type 
used in the glide bomb here described, part of' the values in the coeffi- 
cients a to e (p. 16) have been made practically constant. These 
values will probably not essentially deviate for automatic pilots of 
slightly different form; the larger parts of the values, however, may 
be varied or are known only rather uncertainly so that their influence 
must be taken into consideration. If, moreover, al 1 possible combina- 
tions should be formed, such an enormous number of cases to be calculated 
would result that one will try to attain the goal by a less tiresome 
method. 

Starting from a principal case, characterized on the whole by mean 
values of the variables concerned, the influence of these quantities on 
the lateral stability was investigated one at a time. One forms — if 
the expression be permitted — the partial differential quotients of the 
influence on the lateral stability with respect to all possible vari- 
ations of the controlled automatic pilot. In this way one obtains a 
complete picture of in what respects alterations on the glide bomb must 
be made in order to improve its lateral stability. Of course, there is 
no guarantee that for a different selection of the principal case the 
influences would have made themselves felt in the same way. It is even 
feasible that an influence, favorable for a certain case, might turn 
into the opposite for another. However, where that is the case, one 
will alwayB spot-check anyway, in order to avoid the minor inaccuracies 
of this method. At any rate, the expenditure in calculation for this 
method is the minimum imaginable for a clarification of this rather 
involved problem (with the objective set up initially in mind); thus 
this method was chosen. 
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1. The Stability Domain 

For visualization of the numerical results to be discussed later, 
a form of representation was selected which is also customary in the 
treatment of the lateral stability of uncontrolled airplanes. There 
in a lp,np-plane the domains are plotted where, with the remaining 
airplane values otherwise retained, lateral stability or instability 
prevails. These domains are bounded by curves, the course of which 
shifts if one of these airplane values is altered and appears then as 
a parameter of these curves. 

The value Ip (rolling moment due to sideslip) which is plotted 
on one axis and has for the uncontrolled airplane the significance of 
an indirect restoring force about the longitudinal axis, is like np 
(the yawing moment due to sideslip, or weathercock stability) in 
general controllable within wide limits, and is particularly important 
for the lateral stability. For the controlled glide bomb of this 
investigation lp is replaced by 1 ^ (aileron moment) which is a 
directly effective restoring moment about the longitudinal axis and 
thus becomes comparable to np. Thus the 1_, no-plane logically 
suggests itself for representation of the regions where stability or 
instability of the lateral motion prevails. 

In order to determine whether at any point of the l_,ng- plane 
stability prevails and, therefore, the disturbance of the lateral 
motion shows damping, one must investigate whether all roots X of the 
frequency equation have a negative real part, and whether consequently 
the frequency equation is a Hurwitz equation. This is the case when 
the 5 Hurwitz determinants 

Dj — a 
Dg = ab — c 

*D^ = (ab — c)c — (ad — e)a 

= (ab — c)(cd — be) — (ad — e ) 2 



formed from the coefficients of the frequency equation of the fifth 
degree 

X 5 + ax^ + bX^ + cX^ + dx + e = 0 


are al i positive. 

- -y— — ~r ■ - ^ ' r—r ~ 

NACA Reviewer’s note: D 3 as presented is as in the German text. 

This reviewer believes, however, that in this determinant e = 0 and 
should be used as such. 
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Inversely, if all roots X of the frequency equation hare a 
negative real part, all values Dj_ to D5 are positive. It is a 
consequence of this criterion (and .completely equivalent to it) that 
all coefficients a to e, furthermore D2 and D^, turn out positive. 

If one now visualizes the point in the l^n^— plane considered 

before (which assume dly has been found to be stable) as moving, the 
real parts of the roots, the coefficients a to e, and the 5 Hurwitz 
determinants change. If* one now comes to a point where, for instance, 
the real part of a real root — and thus this root itself — vanishes, 
e too vanishes, since |e| is proportional to the magnitudes of all 
roots j on the other hand, if one reaches a point where the real part 
of a complex root is zero, becomes zero, as can easily be shown. 

In any case, however, D^ = eD^ vanishes, and since for reasons of 
.continuity the totality of all points in the l^n^-plane where stability 
prevails must be connected, the following important statement is valid: 

The limit of the stability domain lies at D^ = 0. 

Thus it is sufficient for graphic representation of the stability 
domain to plot = 0, or, simpler, e = 0 and Djj. = 0 as limiting 
lines of a region at an arbitrary point of which stability is known to 
prevail, as illustrated in figure 4. 


c. Frequencies and Dampings 

In addition to the deter min ation of lateral stability, it is often 
important to know what frequencies and dampings occur when a disturbance 
is damped. Thus one investigates the magnitude of the roots X of the 
frequency equation for a certain point of the l^n^— plane. 

At leant one root is always real and may be numerically determined 
without difficulties. In order to procure an approximate value for a 
root, one considers the problem itself and attempts to make the essen- 
tials of the prescribed motion after a disturbance stand out under 
simplifying assumptions, and thus to obtain an indication for an 
approximate solution. 

According to p. 4, the turning of the flight path into its 
undisturbed direction takes place as an aperiodic transitory motion 
(yawing motion) (fig. 5) on which damped rolling and yawing oscillations 
are superimposed. 

Of course, this applies only under the condition that the restoring 
moments about the longitudinal and vertical axes are so large that 
really oscillations, not aperiodic motions, originate. Then it will be 
possible to replace their influence on the yawing motion by the 
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influence of a mean position of the glide "bomb about which the rolling 
and yawing oscill a tions are performed. By the isolation of these two 
types of oscillations one arrives, therefore, at a consideration of the 
yawing motion alone; for obtai n ing the approximate solution one has only 
to express the dynamic relations corresponding to this train of thought. * 

A motion about a position of equilibrium vanishes — as can easily 
be seen from the differential equation of a free oscillation — if the 
moment of mass inertia and' the damping moment vanish, or if the 
restoring moment by far predominates over all other moments. Under the 
first assumption and further neglect of the coupling moments stemming 
from the velocities of rotation the three equations of motion from 
pp. 11 and 12 then take the following form: 

c a | v^EV + c^* | vSp = - v (cos 7)^ 

L p’P + V? + L k' k = L P* P + L l' £ = 0 

Np'P + + N k «k = + Ng*| = 0 

In the second and third equation the partial moments, separately 
dependent on qp and k were again conqprised into a single moment 

caused by I alone. In this manner it may easily be seen that both 

equations can be satisfied only when p as well as | is permanently 
zero. 


Of course one obtains the same result if one assumes, according to 
the other way of thinking, that the restoring moments by far predominate 
over all other moments since in that case even a vanishingly small 
angle of sideslip P or aileron deflection g produces equilibrium 
of moments . 

Thus with P = 0 there becomes, according to p. 12, cp = \i and 
k = X cos 7 , and according to pp. 7 and 8 , the aileron deflec- 
tion g = k( k + mcp ) = k( X cos y + up) =0 or j a = — c0 ^ 7 x. 

This signifies that during the yawing motion the mean bank is 
proportional to the error in the course of the flight path and is corre- 
spondingly neutralized when approaching the direction of the undisturbed 
flight path (as it was described at the beginning, p. 4, and now con- 
firmed, somewhat mere accurately, by calculation). 
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If this value for p is inserted, into the first equation, there 
results a solution of the remaining differential equation of the first 
order : 

X = -^- 

mTj, 

This root describing the yawing motion is — as was to he expected — 
independent of 1® and np since the magnitude of the frequency of 
the rolling and yawing oscillation about its mean position is, at 
first, of minor Importance. Only in the proximity of the ly 
and np— axis wnere (as was said before) the restoring moments are so 
small that the presuppositions of this consideration are no longer 
satisfied, larger deviations from the otherwise rather accurate approxi- 
mation for X occur, as figure 7 shows for a later example. 

Analogously, one may obtain indications for the frequency and 
damping of the remaining rolling and yawing oscillation. The rolling 
oscillation stands out clearly if one visualizes the disturbances about 
the vertical axis as completely damped; of course, the influence of the 
coupling moments also must be neglected. With 3 = X = 0 there remains 
in the equation ( 2 a), p. 15 , only the first line 


ds 


T- 


F 


V “ ^ ~ “ 0 


From this ordinary oscillation equation then result the approxi- 
mations (required for the later numerical evaluation) for the product II 
and the sum Z of the (mostly) conjugate— complex pair of roots for the 
rolling oscillation; 


Therewith the frequency f x and the damping of the rolling oscil- 
lation become 

■ ih ^-(i ) 2 cps ani - -^(t) 

Therein the damping is explained as the ratio . of the damping factor 
of an oscillating configuration to its critical damping factor; thus it 
is zero for an undamped oscillation, and 1 for the aperiodic limiting 
case. 
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The way to an estimate of the yawing oscillation is perfectly 
analogous. If the disturbances p and x are now put equal to zero, 
one obtains from ( 3 a), p. 15 



and the product ir and the sum a of the pair of roots for the yawing 
oscillation are approximated: 



the frequency f z and the damping D z of the yawing oscillation then 
are 

- 6^3 l/* "(I) 2 CPB - - ¥$) 

Later results (figs. 9 and 11) show that the values II and it 
which are decisive for the frequency are in good agreement with the 
actual values, whereas the values determining the damping £ and tr 
frequently agree with them less the more one approaches the limits 
of the stability domain. 


YI. NUMERICAL EXAMPLES 
a. Numerical. Data 


As mentioned at the beginning, the object of this investigation 
is to give an insight into the lateral stability of the glide bomb 
treated here in order to obtain indications for a design as favorable 
as possible. The required numerical treatment is outlined for this 
special case only; however, even for attempting a quite general solu- 
tion of this problem there is no other way than to use average values. 
The seeming limitation of general validity Is a necessity due to the 
nature of the matter, because its inter-relations are so complicated 
and not readily seen. Nevertheless we shall attempt, by means of 
consideration of relatively extensive variations for the most essential 
properties of this glide bomb, to include all similar automatic pilots. 
However, be it here stressed once more that a comparison is possible 
only when the same method of control is taken as a basis. 

The numerical data used stem partly from wind tunnel measurements 
and partly from rough calculations . Since the automatically-controlled 
bomb no longer has a form remotely similar to an airplane, several 
quantities deviate rather far from conventional ones as will be seen 
from the following conciliation. Therein the values which form the basis 
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for the principal case described on p.19 are underlined; moreover, the 
limiting values within which the stability domains were plotted are 
given. 


Compilation of numerical data 

Lift: c a = 0.15, 0.35# 0.55; in reference to an increase in weight 

(treated separately later), assuming constant velocity, further 
increases in lift values up to 0.643 may result. 


Lateral force: the slope (partly dependent on the variations in magni— 

8c n 

tude of the front stabilizer and vertical fin) is c^* = = 0.3 j 


0.5, 0.7; the Influence of the velocities of rotation a^. and a> z and 
of - the aileron deflection £ were neglected, as mentioned before on 
pp. 9 and 10. 


The following moment coefficients oj 
obtained under the assumption that 


and 


at- 


1.5 to 2.5. 


the lateral stability were 
= 15 to 20 


Dancing in roll: 1^. = 4, T, 10; since it is a function of the fairly 

constant c a », it can be arbitrarily varied (at the most, insignifi- 
cantly, by front stabilizer and vertical fin); nevertheless the influence 
of a possible inaccurate estimate or measurement will be clarified by 
talcing the two limiting values into consideration, too. 


Yawing— rolling moment : l z = 1-45* 2. 45, 3-45; it varies with c a and 

is, moreover, for the principal case (c a = O.35) altered according to 
the given limiting values. 

Rolling moment due to sideslip: lg = 0, 2, 4j aside from the less 

important influence of . the front stabilizer, it may be decreased very 
considerably (according to oral information by A. Schteferdecker) by 
end plates on the wings so that the noteworthy case of a vanishing 
rolling moment due to sideslip is investigated as well. 

Aileron rolling moment: If = -2.6; it is given by the ailerons and has 

for the controlled glide bomb not such an Independent significance as the 
aileron moment: lq, = -kml$ = 0 . . .3; this change for lqj signifies 
for an average value of 0.2 for the measuring ratio m a variation 
for the control— gearing ratio k from 0 to over 5« 
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Bolling— yawing moment : 


it varies with c t > = ^-4, thus is a function 


of a(c a ) and has for c a = O.35 the approximate value n^- = -0.l4; 
since it appears only in the combination (n x — T v y), the Influence of 
the gyroscopic term may be simultaneously taken into consideration in 
this expression by suitable selection of the limiting values 
for (r^ - T F y) = -0.279, -0.140, + O.329. 


Dancing— in— yaw : n z = 2, 4, 6 ; the relatively small damping— in— yaw is 
improved by the front stabilizer and vertical fin and is decreased by 
negative c^, and thus for the larger c a -values (which has to be noted 
when the case arises). 


Yawing moment due to sideslip (weathercock stability): np = 0 . .0.6; 

it is determined by the position of the center of gravity; since the 
latter is mostly fixed, for reasons of longitudinal stability, this 
moment may be altered appropriately by front stabilizer and vertical fin. 


Aileron yawing moment : its ratio to the aileron rolling moment 

is nj: /l^ = —0.01, 0, +0.01; this sometimes uncertain and fluctuating 
value may become inportant for the sign of e (as discussed on p. l8), 
particularly for a larger rolling moment due to sideslip, and is there- 
fore taken into account in spite of its apparent insignificance. 

An increase in the weight of the glide bomb for constant radii of 
inertia and an increase in flight altitude enter into the relative 

aircraft mass density: = -SL. = 700, 1000, 1300; the first value 

gpsF 

corresponds, for an altitude of H = 1.0 kilometer, to a weight of 
of G = 600 kg; the two other values correspond, for this weight, to 
flight altitudes of 4.5 and 6.9 kilometers (calculated according to the 
very accurate rule of thumb by Kholler [lO]) : 

p = 5 ~ ~ jj; ,kgs 2 /n£ 
o 20 + H 

Aerodynamic time unit : Tj. = p s |r = 4.48s for the principal case 

at v = 125 meters per second; it is, at first, not freely selectable but 

is a function of several values already mentioned and of the longitu dinal 

motion of the glide bomb, and thus changes for every case individually. 

The same applies to the longitudinal inclination of the flight 

path: tan 7 = - 0.147 and to the variation with time of the longitudinal 

inclination: y = o7L. — 

■> s 

Eegarding the control one may alter, aside from the control-gearing 
ratio k discussed for the aileron moment lm, the measuring 
ratio: m= 0.1, 0.2, 0.3; for the control used here it depends on the 

angle formed by its measuring axis and the vertical axis of the glide 
bomb which is in these three cases approximately 6°, 11 °, and 17°. 
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b. The Principal Case 

Although the principal case vhich is characterized by a steady 
rectilinear flight path probably never actually occurs, it has its 
Justified significance for an investigation like the present one, 
because it avoids as far as possible any peculiarity in any respect. 

Thus one does not take a certain, actually occurring phugoid motion 
as basis for it; one starts, on the contrary, from the conventional 
concept that this motion has, inversely, to be regarded as deviation 
from a steady rectilinear flight path and has, accordingly, only the 
significance of a special case. 

The principal case, selected according to this point of view, may 
now be treated numeric ally , using the numerical data of the previous 
paragraph. One obtains a good survey always by plotting first the 
stable domain for the representation of which, according to p. 21 f 
the line s e = 0 and 1% = 0 are drawn into a l^np— plane . e = 0 is 
a hyperbola degenerated into a pair of straight lines and can easily 
be plotted into the l^n^— plane . In contrast, 

Djj. = (ab — c)(cd — be) — (ad — e) 2 = 0 is of a higher degree in 2^ 
and n B , and requires, for the numerical evaluation, some deliberation 
in order to obtain with a minimum of calculation expenditure sufficiently 
accurate results. The numerical treatment of such problems' mostly is 
discussed rather cursorily; however, the attainment of fi n al results is 
still a long way off, and thus a few useful remarks concerning the 
calculation will be inserted here. 

First, it is always useful to clarify the desired accuracy in order 
to determine accordingly the number of digits required for calculation. 
Here the result ought to show three digits, the last of which may be 
uncertain; since, however, in the calculation process the first digit 
of many an important number is lost in the forming of differences, four 
digits are necessary to start with. Thus, the rounding up or off did 
extend not only to whole but also to half units of the last digit 
carried which was then denoted by Burrau’s point 11 . This simple 
m Finns permits, by the way, for the same number of digits an effortless 
doubling of the accuracy of a calculation if the necessity arises. Any 
carry ing of further digits would, after all, be wasted effort if they 
finally are not stressed in the graphical representation and do not 
affect the result in any other way. By making suitable proofs during 
the le ng thy calculation process, the always occurri n g un av oidable 
errors may be prevented from doing more extensive damage. 

Figure 6 shows the thus obtained result, the stability domain of 
the principal case. The continued necessary determination of the 
stability at an arbitrary point of the represented I—no— domain outside 
of the negative (shaded) region of e and is quite sirtple: the 

coefficients a to d there directly turn out positive. Just as the 
Eurwitz determinant Dg = ab — c which was already formed in 
calculating D^. 
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Consideration of the principal case as well as of later cases 
shows quite uniformly that the moment coefficients 1^ for the aileron 
moment and np for the yawing m o m e nt due to sideslip (weathercock 
stability) are more or less equivalent with respect to lateral 
stability; in a representation at the same scale for I-) and for n B 
the approximate mirror image symmetry of the regions with respect to P 
the bisector of the angle would become even clearer. A co mp ariso n of 
the results shows that this mirror image symmetry is more or less 
strongly influenced by various couplings and dampings about the longi- 
tudinal and vertical axis. Hence the equivalence of an immediately 
effective restoring moment for these two axes — which had been assumed 
previously in selecting the lmnp-plane — proves correct; under similar 
circumstances it probably holds true for other methods of control as 
well. 


The large stable domain of the principal case can be utilized 
practically only with a few restrictions. It will be useful — if only 
because of the uncertainty (discussed on p. n) in actually maintaining 
a certain value of the weathercock stability — to select only operating 
points sufficiently far distant from the lqj-axis, even if theoretically 
stability prevails up to their immediate proximity. For the aileron 
moment, conditions are not so sensitive; however, one will also select 
values at least ample enough that (in spite of the always possible 
inaccuracies of construction) no excessively erroneous positions remain 
permanently in existence. An upper limit for l<p and np will be 
given mostly by the frequencies which the glide bomb is not to exceed 
with reference to its longitudinal and vertical axis. 

A possible range chosen according to these deliberations might for 
instance lie at 1^ = 1 and n B = 0.3; the stability for this case 
will now be investigated more closely. In order to be able to Judge 
hew it must be shifted if its behavior concerning the initial disturb- 
ance causing the oscillation is to be altered, and in order to gain a 
more exact insight into the distribution of frequencies and dampings in 
the stability domain, the latter are additionally determined along 
two sections at 1^=1 and np = 0.3* 

First the real root of the frequency equation describing the 
yawing motion is determined for the entire l^ynp— region represented. 

As an approximation the value given on p. 23 





0.35 

0.2 x 4.48 



was used, and the further correction was made according to Hewton with 
the aid of Homer’s scheme. 
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Figure 7 shows the result. As was to "be expected according to 
previous deliberations, major deviations of the root X from the above 
approximate value, otherwise accurate within 1 to 2 percent, occur in 
the prox imi ty of the 1 — and n^— axis . This approximate value signifies 

that accor ding to = 1.8s an error in course has been damped to 

“A 

half amplitude. 


If one can make sure, in stability investigations of this or a 
similar kind, that one stays within the part of the lg^ip— plane only 
where this good agreement exists, it is very profitable beforehand to 
lower the degree of the frequency equation by one with the aid of this 
root. As a result, one need calculate only with much simpler Eurwitz 
determinants and the evaluation of the frequency equation is facilitated. 
Unfortunately, this method Is not possible for the present investigation, 
because the stability domains are, for the sake of a complete survey, 
plotted to the 1 — and n^— axis. 


If the root of the yawing motion has been determined with suffi- 
cient accuracy for an arbitrary point of the l_np— ; plane and is inserted 
for the last time as a proof into Homer's scheme, there result at the 
same time the coefficients a' to d' of the remaining frequency 
equation of the fourth degree 

X + a'x3 + b'X 2 + c'X + d* =0 


Its roots describe the oscillations of the glide bomb about the 
longitudinal and vertical axis which in the proximity of the Imp- 
end n p-axi s always become aperiodic; this equation now has to be solved 
numeric ally . The iteration method of v. Karman— Trefftz [4] has the 
disadvantage that it converges poorly for roots of equal order of magni- 
tude an d for that reason is not always usable. In order to attain the 
aim quickly also for this case, an approximation method is given in the 
appendix which permits a rapid and clear solution. Figures 8 to 11 
were calculated according to this method — chiefly according to the 
"pi — method" because It agrees better with the approximate values 
(pp. 23 and 24) used generally. 

Figure 8 shows the section along 1=1 through the stability 
domain of the principal case. No peculafities of any kind occur. The 
rolling oscillation remains almost unchanged; the yawing oscillation 
becomes slower and slower for smaller n^— values and Its two roots 
finally become real so that it Is no longer noteworthy in this boundary 
region for vanishing yawing moment due to sideslip at the limit of 
static stability. The damping of the yawing oscillation meanwhile 
increases more and more rapidly and would exceed the value 1 in the 
aperiodic limiting case. 
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The approximate values for U and n agree well, those for z 
and a only Just tolerably, with the exact ones, as can he seen from, 
figure 9* 

The section along np =0.3 through the stability domain of the 
principal case shows a somewhat different pioture insofar as one 
approaches here, with diminishing aileron moment, the limit of dynamic 
stability, as can also be seen from figure 6. Thus one damping vanishes 
there, in this case that of the rolling oscillation. No further points 
were calculated below the value 1® — 0.3, where this occurs, since they 
have no practical significance. The fact that the damping, which finally 
disappears at the stability limit, starts decreasing more or less 
strongly before that, offers another clue for the selection of the 
operating point where the dampings — not too good as it is — are to be 
utilized as efficiently as possible. 

The agreement of the approximations with the exact values of IT, 

Z, it, and a can be seen from figure 11, which again shows the good 
serviceability of the approximate values of IT and it determining 
the frequencies, and the not so good serviceability of those of Z 
and cr decisive for the dampings. 

What was said above concerning the root X of the yawing motion 
and the variation of the frequencies and dancings applies, with the 
required changes, to any other of the cases discussed below that result 
from the principal case, and need not be repeated. 


c. Influences of the Glide Bomb Values 

Starting from the principal case one now varies in turn one 
quantity, as described on p. 19, and investigates its influence on 
the stability domain. All values used in this problem were thus 
included in arbitrary sequence, no matter whether their variation can 
be achieved easily or not at all; reasons for this will be given later. 

First, the influence on the stability domain is investigated for 
steady rectilinear flight performed with different c a and accordingly 
varied velocity, the other ratios remaining unchanged. As always in 
uncontrolled lateral motion here also a distinct c a dependence becomes 
apparent; for smaller c a — values the stability domain increases. 

For the important determination of the variation of the stability 
domain during a-phugoid oscillation one made the assumption that the 
release velocity of the glide bomb is only 60 percent of its equili- 
brium velocity; thus it was based on a very pronounced path oscillation 
which shows at the four locations denoted in figure 13 the following 
roughly calculated value; 
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TABLE 1.- VALUES LURING THE PHUGOID 


Location 

B 

7 

. l 

(n^. - Tg-r) 

1 

wave crest 

75 

-0.146 

-O.062Q 

O.329 

2 

descending flight 

125 

-. 589 

0 

0 

3 

• 

1 

3 

wave trough 

156 

-. 146 

.0386 

-.279 

4 

ascending flight 

125 . 

.126 

0 

-.140 

steady rectilinear flight 

125 

-0.146 

0 



-O.i4o 


Flight in the wave trough is found to he most unfavorable for later 
lateral stability. The considerable difference between the stability 
domains of the wave crest and the wave trough stems solely from the 
additional gyroscopic moment due to the rotation of the path of the 
longitudinal motion and shows clearly that the taking into consideration 
of the gyroscopic term which was motivated on p. 6 is justified. Since 
it appears only together with the rolling-yawing moment n x in the 
combination (n^. — T -$ 7 ), — as both vary with — the two stability 
domains for the wave crest and the wave trough may be interpreted also 
as if in the steady rectilinear flight of the principal case the n^. 
were one time equal to +O.329 and the other time equal to —0.279* Thus 
the influence of the rolling— yawing moment on the stability domain is 
simultaneously clarified and need not be investigated separately. 

Compared to the locations 1 and 3 where the path rotation of the 
longitudinal motion is most pronounced, the approximately linear parts 
of the phugoid motion at the locations 2 and 4 exert hardly any influ- 
ence on the stability domain. The longitudinal inclination of the 
flight path 7 as such, is, therefore, insignificant provided the 
c a — value remains the same. 

An actually performed flight of a glide bomb with the approximate 
values of the principal case is represented in figure 14. The aileron 
moment is about lq> = 1, the weathercock stability np = 0.5; it is 
true, the reduction in velocity at the release is not as high as had 
been presupposed in figure 13 so that the phugoid motion appears less 
marked. The lateral curvatures of the path at about 30 to 40 and 50 
to 60 seconds flight duration stem from arbitrary commands which are to 
be interpreted as disturbances for the lateral motion in the sense of 
p. 5* 

















32 


H ACA TM 1248 


Another value 'besides the control— gearing ratio k which is vari- 
able by the control is the measuring ratio m; for the type of control 
used here it depends on the tilt of its measuring axis in the glide 
bomb. Figure 15 shows that in the case m = 0.1 (too slight tilt) a 
large part of the useful stability range is lost because one is already 
too close to the limiting case m = 0 where a bank is no longer 
measured and the control method is therefore useless. Also, for m = 0.3 
(excessive tilt) conditions are somewhat more unfavorable than 
for m = 0.2. Thus an optimum regarding the magnitude of the stability 
range lies between m = 0.1 and m = 0.3, probably in the proximity of 
the value used for the principal case m = 0.2. Since the time to damp 
half the anplitude of the yawing motion varies with m, as can be seen 
from pp. 28 and 29, it is from this point of view expedient to use a 
small m, thus a small tilt. 

The influence of the rolling moment due to sideslip may be best 
understood on the basis of the conventional representation of the 
stability range for the uncontrolled lateral motion in the lpnp— plane 
as it can be found for instance in Mathias [12] . For this purpose the 
representation of the stability range for the principal case without 
control, thus for lqj = 0, was incorporated in figure 16. If one 
considers in it, for instance, the point 0.3 on the no-axis, one is in 
the statically unstable range and, as is well known, it takes a definite 
minimum rolling moment due to sideslip Ip to attain stability. If 
one proceeds in the direction of increasing lp— values, one reaches — it 
is true, for rolling moments due to sideslip so large that one would 
no longer consider them, for instance, for airplanes of ordinary 
design — again an unstable region, which this time, however, is the 
region of dynamic instability. If one stops at lp = 4, one has reached 
approximately the rolling moment due to sideslip of the automatic pilot 
considered here, which is, therefore, much too large. In order to 
attain stability, the control must intervene with a minimum aileron 
moment I<p, which would be read off for np =0.3 from figure l6a. 

For lp = 2 and the same np the glide bomb with the values taken here 
as a basis would be stable also without* control, and only below np = 0. 15 
a control would become necessary. The case is slightly different 
for lp = 0. Here the uncontrolled glide bomb is not stable for any 
positive np; however, a very small 1^ on the part of the control 
system is sufficient to eliminate this statio instability; hence in 
figure l6a the different character of the unstable range for lp = 0 
and the two other lp— values. It is shown with particular enphasis how 
many advantages a small rolling moment due to sideslip offers, which 
probably applies likewise to glide bombs of similar type. 

The weight of the glide bomb may be increased for unchanged lift 
coefficient c a and correspondingly increased flight velocity or 
inversely. The first case is represented under the assumption that the 
radii of inertia remain unchanged; it shows that by using this method 
surprisingly little of the stability range is lost. 
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As discussed an p. 26, the same over-all picture would hold if for 
unchanged weight the flight altitudes were instead of the original 
110 kilometer (principal case) now 4.5 and 6.9 kilometers. This fact 
signifies only that the stability range remains fairly constant when 
large differences in altitude are passed through in flight. 

Figure 18 represents the flight of a glide bomb of about 1000 kilo- 
grams weight to which corresponds, for an altitude of 2.0 kilometers, 
a relative aircraft mass density n s = 1300. The initial lateral curva- 
ture of the path was not intentional but was probably due to release 
disturbances. For a prescribed size of the glide bomb, the high flight 
velocities of more than 200 meters per second occurring here can be 
avoided only by an increase of the lift coefficient c a , as had been 
presupposed for figure 19* 

Here an increase in weight, for unchanged radii of inertia, is 
obtained at the price of a much more extended enlargement of the 
unstable region, as had been the case for unchanged lift coefficient c a . 
Comparison with figure 12 where only the variation of c a was investi- 
gated, while the weight remained unchanged, shows a striking similarity 
in the unfavorable influence of high c a — values on the stability region. 
Thus the assumption suggests itself that of the possible combinations 
of lift coefficient, flying weight, and velocity the first exerts the 
decisive influence on the size of the stability domain. In order to 
recognize a possibly existing connection, the table 2 was set up: 


TABLE 2.- RATIO OF THE SIZE CF THE UNSTABLE REGION AND 
DIFFERENT STATES OF FLIGHT FOR RECTILINEAR PATH 


Figure 

Size of the 
unstable region 

c a 

Us 

Fr 

T f s 


0.014 

0.15 

700 

2340 

2.93 

12 

.117 

• 35 

700 

1000 

4.48 


.270 

• 55 

700 

610 

5.66 

23 location 2 
location 4 

0.122 

0.35 

700 

1000 

4.48 

.120 

• 35 

700 

1000 

4.48 

IT 

0.142 

0.35 

1000 

1430 

5*35 

.162 

• 35 

1300 

1850 

6.11 

19 

0.259 

0.497 

1000 

1000 

6.4o 

• 397 

.643 

1300 

1000 

8.32 
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The rectilinear but not steady parts of the ascending and descending 
flights during the phugoid motion also were included with the cases suit- 
able for this coup arisen. As measure of comparison for the unstable 
region we selected the size of its area in the represented lmnp— region, 
expressed in ( l<pnp )— units . The weight was expressed by the relative 
aircraft mass density u s and the flight velocity by the Froude 
number Fr = v^/bg. Finally, the aerodynamic time unit Tp = p s ^ 
was included sinoe it belongs in this example. 

As even a most superficial consideration reveals, the magnitude of 
the unstable region is an approximately unequivocal function only of the 
lift coefficient c a j this dependence is represented in figure 20. 

Without generalizing too hastily it is, therefore, shown with 
satisfactory regularity that with respect to the magnitude of the 
unstable region, small o a — values are favorable, large ones unfavorable. 
Accordingly, it will be best, if a minimum unstable region is desired, 
to fly with a small c a — value, mostly given by the maximum admissible 
flight velocity. If the longitudinal motion of the glide bomb is 
controlled arbitrarily or automatical 1 y, the unfavorable influence of 
large angles of attack maintained for a longer time must be taken into 
consideration in a given case. 


At this point it should be mentioned explicitly that in oase of 
a change of c a the two related coupling moments l z and n x were 
also changed throughout, because the same holds true actually. The 
Influence of the rolling moment due to yawing l z as such is insignif- 
icant, as figure 25 will show, and the rolling moment due to yaw n^ 
also is only of secondary importance compared to Oa- It must be borne 
in mind that in figure 13 an n x that should appertain to c a = 0.70 
corresponds to the slightly enlarged unstable region of the wave trough. 
Hence it results that it is really the lift coefficient c a , and not one 
of the coupling moments l z and in. dependent on it, which exerts the 
decisive influence on the size of tne unstable region. 


Figure 2la shows the expected favorable influence of a large 
dartping— in— roll. For 1^. = 4, where the moment coefficients of the 
damping— in— roll and damping— in— yaw are equal, the mirror image symmetry 
with respect to the bisector of the angle mentioned on p. 28 is strikingly 
good,* this stands out clearly in figure 21b where l<p and np are 
plotted to the same scale. , If the damp ing-ii>-r oil is improved, the 
unstable region is reduced in such a manner that it decreases more in 
dtp— than in np— direction (again under the assumption of equal scale! ) 

The variation of the damping— in-yaw now also has the expected 
Influence, as will be. shown by the following oase. 
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Here, too, the self— evident favorable effect of large damping— in- 
yaw is obvious j the unstable region is influenced mainly in its 
extension parallel to the np-axis. 

One is not always able to make the damping— in— yaw — and even less 
the dashing— in— roll! — so large as to obtain a decisive effect on 
lateral stability. The investigation of these and also of other cases 
where an influencing of the aerodynamic properties of the glide bomb is, 
in practice, hardly possible has no other purpose than to determine 
which one of the two limiting cases of numerical data known only 
inaccurately must be substituted into the calculation to stay on the 
safe side. 

The lateral force which, among other effects, also has a damping 
effect, exerts a surprisingly small influence, the reason being that it 
is by one order of magnitude smaller than the moment coefficients of 
the damping— in— roll and damping-in— yaw. 

For negative n| yl^ , that is, when the restoring force about the 
vertical axis resulting from the weathercock stability is reduced 
because of the aileron deflections simultaneously produced by the 
control, the coefficient e may become negative in spite of positive 
(but too small) weathercock stability. This case discussed on 
pp. 17 and 18 is practical ly the only one where static lateral insta- 
bility also is possible. 

For positive n| /I4 , on the other hand, the region of dynamic 
instability is enlarged so that — the aileron yawing moment being not 
definitely known — it is best to take both influences into consideration. 

The rolling moment due to yawing has praotloa.il, y no influence at 
al 1 on the magnitude of the stability domain, as is to be expected 
(if for no other reason) because of the very full wing shape of the 
glide bomb. It is related to co z , exactly like the gyroscopic moment 
that would arise for a path rotation of the longitudinal motion and, 
accordingly, would be just as insignificant. Thus it is shown also 
from this point of view that the neglect of the gyroscopic term in the 
equation of moments for the longitudinal axis was fully justified. 


d. Conclusion 

In surveying the figures which represent the stability . domains 
as functions of the various glide bomb values, one finds that 
directio nal stability is with certainty attainable for every case, by 
appropriate adjustment of the control in cooperation with the 
ailerons (m,l<p) and a corresponding yawing moment due to si-deslip (np) . 
No other definite aerodynamic properties of the system need be assumed 
to make the control method treated here serviceable although some are 
more or less unfavorable. Of course, one will, where it is easily 
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possible, select tile more advantageous design; the following table may 
be used as a directive. 


TABLE 3.- INFLUENCES CUT THE LATERAL STABILITY 


Parameter 

Figure 

Favorable 

Unfavorable 

Lift coefficient c 

a 

12 

Small 

Large 

Phugoid 

13 

Wave crest 

Wave trough 

Rolling— yawing moment n x 

13 

Large ( > 0) 

Small ( < 0) 

Measuring ratio m 

(a function of gyroscope— 




tilt angle) 

15 

About 0.2 

Other values 

Rolling moment due to 

16 



sideslip lp 

Small 

Large 

Weight increase (p s ) 

17, 19 

Small c a 

Large c a 

Altitude H 

17 

Small 

Large 

Damping— in— roll 1^ 

21 

Large 

Small 

Davping— In— yaw n z 

22 

Large 

Small 

Lateral— force increase c^ 

23 

Large 

Small 

Aileron yawing moment ng 
Rolling moment due to yawing l z 

24 

About 0 

Other values 

25 

Practically without 
influence 


The table is compiled without consideration of the question whether 
or not an alteration is actually feasible for the individual case; at 
the same time, the table should be used to select, in case of unreliably 
known values, always the most unfavorable values for safety reasons. 

Most significant with regard to the minimum possible unstable region 
for lateral motion of the investigated controlled glide bomb are small 
rolling moment due to sideslip lp and small lift coefficient c a . 

As a conclusion it should be mentioned again, in regard to what 
was said on p. 5 , that the object of this investigation was not to 
present, for a certain Individual case, a calculation as accurate and 
complete as possible, but rather to give a survey so comprehensive that 
Its results make a more general insight into the problems of this 
automatically controlled glide bomb possible. 
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YII. SUMMARY 


The investigation of tie lateral stability of an automatically 
controlled glide "bomb led also to the attempt of clarifying the influence 
of a phugoid oscillation — or of any general longitudinal oscillation — 
on the lateral stability of a glide bomb. Under the assumption that 
its period of oscillation considerably exceeds the rolling and yawing 
oscillation and that c a is, at least in sections, practically constant, 
the result of this test is quite single. It becomes clear that the 
influence of the phugoid oscillation may be replaced by suitable vari- 
ation of the rolling— yawing moment on a rectilinear flight path instead 
of the phugoid oscillation. If the flying weight of the glide bomb of 
unchanged dimensions is increased, an increase of the flight velocity 
will be more favorable than an increase of the lift coefficient. The 
arrangement of the control permits lateral stability to be achieved in 
every case; a minimum rolling moment due to sideslip proves of great 
help. 
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YIII. APPENDIX 

AN APPROXIMATION METHOD FOR SOLUTION OF ALGEEGRAIC 
EQUATIONS OF THE FOURTH DEGREE 


One imagines the four roots of a prescribed algebraic equation 
of the fourth degree 

+ a'X 3 + b*X 2 + c»X + d' = 0 

grouped into two pairs; let the root products and sums for each of theBe 
pairs be denoted by IT, E and it, ex, respectively. If two conjugate- 
complex roots exist, their combination is self-evident and real roots 
are treated as pairs, grouped arbitrarily by twos. 

The following ratios are known to exist between the coefficients a' 
to d* and the root products and sums: 


a* = - Z - a 
b* = II + ZG + n 
c« = - Ilex - JtZ 
d« = Hit 


They readily permit the following method of solution for the four new 
real unknowns II, Z, Jt, and a to be read off: 

For instance, let an approximate value of II (or xc) be known; 
then there results first, from the fourth equation, xt (and H, respect- 
ively) and then, as can be seen easily, from the combined first and 
third equation, the cx and Z. If the H (or Jt) used initially 
had been correct, the result (in case of substitution of all 
these successively found values into the right side of the second 
equation) would be exactly b 1 ; but since an approximation had been 
used, the proof II + Iff + it — b* will be omitted. Thus the 

result of the proof is in a very simple manner dependent on H (or xt); 
this value can now be easily determined according to the regula falsi 
(bracketing) with arbitrary accuracy in such a way that the proof gives 
zero; then one has found the correct value of n (or Jt). Thus the 
course of calculation appears as follows, starti n g from an approximate 
value for H: 

II; Jt = — ; cr = — Z = —(a 1 + ct) ; II+Zff+Jt — b 

H II — Jt 
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Every further line is used to obtain a more and more exact value 
of IT until the desired accuracy is reached. If one finally inserts 
the correct value for IT, one obtains, according to this method, not 
only the remaining values it, c, and E, but has in addition for all of 
them a common proof. If one starts with it instead of IT, the calcu- 
lation takes, after determination of IT = d’/*, exactly the same course. 
For IT and it being approximately equal, a more accurate determination 
of the second value of each line is sufficient to have, in forming the 
difference U — it, enough digits left to calculate a. 

It is true that this method will fail in case of IT = it; however, 
the approximation method may be used here also, provided one starts, in 
a slightly varied manner, with E (or cr) instead of H (or it). With 
the aid of the four ratios used initially the following second type of 
solution is Just as clear: 

After selection of E (or a), the a (or £, respectively) immedi- 
ately results from the first equation, and then, from the second and 
third together, the it and IT; the proof Hit — d' ^ 0 now serves 
exactly as in the method described before for Improvement of the initial 
value E (or cr). The separate lines of the calculation then read, if 
one starts, for instance, with E: 


E; cr * -(a* + e); « 


o* + <r(b» — Ed) . 
-(E - cr) J 


H = b» - Ecr - it; Hit - d* 


This "Sigma-method" also has a weak point, in the case of E = cr. 
However, in that case it is very well supplemented by the "Pi-method", 
except for the (practically highly improbable) possibility of simul- 
taneous equivalence II = it and E = o . Anyway , this fact would be 
shown by definite ratios between the coefficients a* to d', and the 
solution is then always possible in some other manner. 

With the values H, E, it, and a the problem is practically 
solved, since the four desired roots now can be easily determined from 
the two quadratic equations 

X 2 - EX, + H = 0 and X 2 - aX + it = 0 


Thus only the character of the roots — whether real or complex — emerges; 
this knowledge is, therefore, not even required for performing this 
approximation method. 


Translated by Mary L. Mahler 
National. Advisory Committee 
for Aeronautics 
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Figure 1.- The subject glide bomb (scale 1:20 dimensions in cm.) 
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Figure 2.- The path oscillation (phugoid) . (Not true to scale.) 



Figure 3.- Concerning the lateral motion of the glide bomb. 
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Figure 10.- Frequencies and dampings (np = 0.3). 
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Figure 11.- TT,£, tt, a and their approximations (n p = 0.3). 
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Figure 12.- Influence of the lift coefficient c on the stability domain. 
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Figure 13.- Influence of the phugoid oscillation on the stability domain. 
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Figure 14.- Path and velocity of a glide bomb with (i s = 700 (principal case). 
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Figure 15.- Influence of the measuring ratio m on the stability domain. 
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Figure 16(a) and (b).- Influence of the rolling moment due to sideslip 1 

on the stability domain. 
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Figure 11 .- Influence of the weight (c & constant) on the stability domain. 
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Path and velocity of a glide bomb with n„ = 1300. 
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Figure 23.- 


Influence of the increase of the lateral force 
the stability domain. 
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Figure 24.- Influence of the aileron yawing moment ng on the stability domain. 



Figure 25.- Influence of the rolling moment due to yawing on the stability domain. 
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